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GASEOUS-HELIUM REQUIREMENTS FOR THE DISCHARGE OF LIQUID HYDROGEN 


FROM A 1.52-METER- (5-FT-) DIAMETER SPHERICAL TANK 
by Robert J. Stochl, Joseph E. Maloy, Phillip A. Masters, and Richard L. DeWitt 

Lewis Research Center 

SUMMARY 

An experimental investigation was conducted to determine the effects of various 
physical parameters on the pressurant gas (gaseous helium) requirements during the 
pressurization and expulsion of liquid hydrogen from a 1.52-meter- (5-ft- ) diameter 
spherical tank. The experimental results were compared with results predicted by, a 
previously developed analytical program which was revised and extended for this inves- 
tigation. Tests were conducted for a range of liquid outflow rates, pressurizing rates, 
and initial ullage volumes at a nominal operating pressure of 34. 47x10^ newtons per 
square meter (50 psia) using nominal inlet gas temperatures of 183, 267, and 322 K 
(329°, 480. 6°, and 580° R). Data were obtained using a hemisphere injector. 

The experimental results show that the inlet gas temperature has a strong influence 
on the actual pressurant gas requirements. There is an average decrease in pressurant 
gas requirements of 1.3 percent for a 10 K (18° R) increase in inlet gas temperature. 
The analytical program was able to predict the actual pressurant requirements to within 
a maximum error of 11.5 percent for all cases. 


INTRODUCTION 

During the past several years, a great deal of effort has been devoted to the prob- 
lems associated with the pressurized discharge of a cryogenic liquid from a tank. The 
main objectives of these efforts have been toward the optimization of a propellant tank 
pressurization system. One phase of this optimization is a precise determination of 
pressurant requirements for any given set of operating parameters (i.e. , tank pressure, 
inlet gas temperature, liquid outflow rate, tank size, etc. ). This knowledge would allow 
the design of a system based on only the weight of gas necessary to accomplish the 



Several investigators have developed analyses (e.g. , refs. 1 and 2) which attempt to 
predict, according to a selected set of simplifying assumptions, the quantity of pressur- 
ant gas required during the pressurized discharge of liquid hydrogen (Lt^). Some of 
these simplifying assumptions may, for certain conditions (i.e. , for various injector 
geometries, tank shapes, and tank sizes), limit the capability of the analysis to predict 
pressurant requirements accurately. Because of these limitations, the validity of the 
analytical results have to be largely based on correlations with experimental results. 

An experimental investigation was started at Lewis Research Center to determine 
the effects of various physical parameters on the pressurant gas requirements during the 
pressurization and expulsion of liquid hydrogen from propellant tanks. The experimental 
results were also used to extend the capability of the analysis of reference 1 to predict 
the pressurant requirements for tanks of general size and geometry under various oper- 
ating conditions. Some results of this program, which used hydrogen gas as the pres- 
surant, are reported in references 3 to 5. In general, the experimental results reported 
in references 3 to 5 indicate that the temperature of the inlet gas, the injector geometry, 
and the tank pressure have a strong influence on the pressurant gas requirements. 

There was between a 1.2 and 1. 4 percent decrease in pressurant requirement for a 10 K 
(18° R) increase in the temperature of the inlet gas depending on liquid outflow rate and 
tank size. These data also indicate that the use of a straight pipe injector could reduce 
the pressurant requirement by 9 to 35 percent, depending on liquid outflow rate and tank 
size, compared with that obtained using diffuser-type injectors. It was also shown that 
pressurant requirements increased in proportion to increased tank pressure level. 

The analysis of reference 1 was used to predict the pressurant requirements in ref- 
erences 3 to 5. It was modified in references 4 and 5 to accommodate a spherical tank 
geometry and also to include the heat transfer from the gas to the liquid surface (which 
was neglected in the original analysis). The analysis predicted the actual pressurant r e- 
quirements for the various operating parameters within 12.4 percent with the exception 
of the tests performed with the straight pipe injector. The analysis overpredicted the 
pressurant requirements for the straight pipe injector by as much as 45 percent. 

If only pressurant gas weight is considered, it appears advantageous to use hydrogen 
(because of its low molecular weight) as the pressurant for the expulsion of liquid hydro- 
gen. Since gaseous hydrogen cannot be used as the pressurant for most oxidizers and a 
dual pressurization system would probably result in an additional payload penalty, it 
would appear desirable to use a single pressurant gas which is compatible with liquid 
hydrogen and any oxidizer. To meet this requirement, helium is the only logical choice 
because it has a low molecular weight and, since it is an inert gas, it is compatible with 
any fuel- oxidizer combination. 

This report presents the results obtained using gaseous helium as the pressurant for 
the expulsion of liquid hydrogen. The use of helium as a pressurant in the expulsion of 


2 



r 


liquid hydrogen adds further complication to both the analytical and experimental pro- 
grams. The ullage volume for this system contains a binary gas, part of which is con- 
densable and part is not. The analytical model must either (1) account for the possibil- 
ity of a binary gas in the ullage volume by diffusion and/or evaporation or (2) assume the 
contribution of hydrogen to the ullage volume to be small enough that the ullage can be 
considered at 100 percent helium at all times, except at the interface where it is con- 
sidered to be 100 percent hydrogen. The analysis used in this report is identical to that 
used in reference 4 with the exception that this analysis treats the ullage volume as if it 
contained pure helium. (For convenience, the modifications to the analysis of ref. 1 
that were used in ref. 4 are presented in appendixes A, B, and C. ) Experimentally 
measured gas concentration gradients in the ullage volume were used to obtain mass and 
energy balances and to verify the 100 percent helium model used in the analysis. 

The primary objective of the test work described herein was to obtain experimental 
values on the helium pressurant gas requirements during the initial pressurization and 
expulsion periods for the various operating parameters and to compare these values with 
predicted ones. A secondary objective was to obtain experimental information on tank 
wall heating, liquid heating, residual ullage energy, and mass transfer in order to gain 
an insight into the reasons for any variations in pressurant gas requirements. 

The tests were conducted using the same two 1. 52-meter- (5-ft- ) diameter spherical 
aluminum tanks that were used in reference 4. The pressurant used for this series of 
tests was gaseous helium. The main test variables were (1) tank wall thickness (0. 762 
and 0. 409 cm; 0.30 and 0. 161 in. ), (2) nominal inlet gas temperatures of 183, 267, and 
322 K (329°, 480. 6°, and 580° R), (3) liquid outflow rates between 0. 282 and 0. 824 kilo- 
gram per second (0„ 623 and 1.817 lb/sec), and (4) initial tank ullages of 4, 28, 55, and 
75 percent. Data were also obtained for pressurization of the tank at various rates 
(3.24x10^ to 11.37x10^ (N/m^)/sec; 0.47 to 1.65 psi/sec) from 1 atmosphere to the de- 
sired operating level. All tests were performed at a nominal operating tank pressure of 
34.47x10^ newtons per square meter (50 psia). 


SYMBOLS 


A area, m^ (ft^) 



C orifice coefficient 

Cjj effective perimeter of interior hardware, m (ft) 

C specific heat at constant pressure, J/(kg)(K) (Btu/ (lb) (°R)) 
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Gr 
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H 
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AM 
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N 
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N l" N z 

Nu 


n or i 
P 


AP 

AP* 


Pr 


specific heat at constant volume, J/(kg)(K) (Btu/(lb)(°R)) 

specific heat of tank wall, J/(kg)(K) (Btu/(lb)(°R)) 

a - aw - T — At -5L 

Vw C w 

orifice diameter, m (ft) 


Z 2 AX P' - P Z 1 AX 
Z At P’ Z At 
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Grashof number, — - P ^ 

M 2 


2 2 

gravity acceleration, m/sec (ft/sec ) 

enthalpy, J (Btu) 

specific enthalpy, J/kg (Btu/lb) 

convective heat-transfer coefficient, J/ (m^)(K) (sec) (Btu/(ft^)(°R)(sec)) 

thermal conductivity, J/(m)(K)(sec) (Btu/(ft)(°R)(sec)) 

flow length, m (ft) 

thickness, m (ft) 

mass, kg (lb) 

differential mass, kg (lb) 

mass flowrate, kg/sec (lb/sec) 

molecular weight, kg/(kg-mole) (lb/ (lb- mole)) 

ideal pressurant requirement, kg (lb) 

number of volume segments 

number of data points used in defining average deviation 

particular volume segments 

Nusselt number, h c L/k 

summing index 

pressure, N/m (lb/in. ) 

2 2 

differential pressure, N/m (lb/in. ) 
orifice AP 

Prandtl number, p C /k 
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Q heat transfer, J (Btu) 

Q heat-transfer rate, J/sec (Btu/sec) 

Q* specific heat- transfer rate, J/(kg)(sec) (Btu/ (lb) (sec)) 

n n 

q heat-transfer rate peir unit area, J/(m )(sec) (Btu/(ft )(sec)) 
R gas constant, J/(K)(k:g-mole) (Btu/ (°R) (lb- mole)) 
r radius, m (ft) 

Ar increment of radius,, m (ft) 

Re Reynolds number, 7LVp/ M 

AS increment of are length, m (ft) 

T temperature, K (°R) 

AT differential temperature, K (°R) 

Tg temperature at the edge of thermal boundary layer, K (°R) 
t time, sec 

At time increment, sec 

U internal energy, J (Btu) 

AU differential energy, J (Btu) 
u specific internal energy, J/kg (Btu/lb) 

V volume, m^ (ft^) 

AV volume increment, m (ft ) 

V velocity, m/sec (ft/sec) 

O O 

v specific volume, m /kg (ft /lb) 

W work, J (Btu) 

X percent of gas by volume 
X .n number of net points in the ullage 
x coordinate in direction of tank axis, m (ft) 

Ax space increment, m (ft) 

V expansion factor 

y thickness within the boundary layer, m (ft) 

Z compressibility factor 

z elevation or vertical distance along tank wall, m (ft) 
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X 

M 

P 


h„ At 

1 + — 2 — 


^w^w^w 



coefficient of thermal expansion, 1/K (1/°R) 
specific heat ratio 

finite increment, or total boundary layer thickness, m (ft) 
latent heat of vaporization, J/kg (Btu/lb) 
viscosity, kg/(m)(hr) (lb/ (ft) (hr)) 
density, kg/m 3 (lb/ft 3 ) 

' M at M»r 2 / C b P 


Subscripts: 

A analytical results 

ad adiabatic 

E experimental results 

f final state or condition 

G gas added to tank 

H internal hardware 

He helium 

Hg hydrogen 

i initial state or condition 

i — f from initial to final state or condition 
L liquid 

n summing index 

o condition prior to ramp 

S liquid surface 

sat saturation 

T total quantity 


6 


t transferred 

tptg times 1 and 2 
U ullage 

w wall 


Superscripts: 

' time index 

* evaluation may be performed at beginning or end of time interval 


Constants for Beattie- Bridgeman equation: 

A (N)(m 4 )/kg 2 (ft 4 /lb) 

a m 3 /kg (ft 3 /lb) 

B m 3 /kg (ft 3 /lb) 

b m 3 /kg (ft 3 /lb) 

C (m 3 )(K 3 )/kg ((ft 3 )(°R 3 )/lb) 

W X H 2 % 2 + X He%e 

X H 2 C H 2 + X He C He 

vT 3 


APPARATUS AND INSTRUMENTATION 
Facility 

All tests were conducted under vacuum inside a 7. 61-meter- (25-ft-) diameter spher- 
ical vacuum chamber (fig. 1) to reduce the external heat leak into the propellant tank to a 

_7 

low value. The vacuum capability of this chamber was approximately 8x10 mm Hg. A 
general schematic of the test tank and associated equipment is shown in figure 2. A heat 
exchanger and blend valve subsystem capable of delivering gaseous helium at any desired 
temperature between 167 and 405 K (301° and 729° R) at a maximum flow rate of 

_o 

9. 98x10 kilogram per second (0.22 lb/sec) was used to control pressurant gas inlet 
temperature. A ramp generator and control valve were used for controlling the initial 
rate of pressurization of the propellant tank. A closed- loop pressure control circuit was 
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Figure 1. - 7.61-Meter- (25-ft-) diameter vacuum chamber. 


used to maintain constant tank pressure during the expulsion period. The liquid outflow 
rate was controlled by remotely operated variable flow valves. The liquid hydrogen out- 
flow from the tank was returned to the storage Dewar. 

Liquid outflow rates were measured using a turbine-type flowmeter located in the 
transfer line. The flowmeter was calibrated (within an estimated uncertainty of ±1/2 per- 
cent) with liquid hydrogen over the expected range of flow rates. The calibration was 
performed at Lewis. Pressurant gas inlet flow rates were determined by the use of an 
orifice located in the pressurant supply line. Tank, line, and differential pressures were 
measured with bonded strain -gage -type transducers (estimated uncertainty, ±1/4 per- 
cent). 


Test Tanks 

The experimental work was conducted in 1.52-meter- (5-ft- ) diameter bare wall 
spherical aluminum tanks. One tank (fig. 3) had an average wall thickness of 0.762 cen- 
timeter (0.30 in. ). The other tank was identical except the wall was chem-milled down 
to an average thickness of 0. 409 centimeter (0. 161 in. )» One stainless steel lid served 
both tanks. The lid housed the inlet and vent pipes and the electrical connections for all 
internal tank instrumentation. The lid was 0. 457 meter (18 in. ) in diameter, 3. 18 centi- 





Gaseous 



valves 


Figure Z - General schematic of facility. 








Figure 3. - 1. 52-Meter- (5-ft-) diameter spherical test tank. Tank wall thickness, 
0.762 centimeter (0.30 in.). 


meters (1. 25 in. ) thick, and weighed 67. 50 kilograms (149 lb). The inner surface of the 
lid conformed to the contour of the tank wall. 

A viewport and television camera were installed on the 0. 762 -centimeter - 
(0. 30-in. -) wall test tank (fig. 4) for the purpose of visually locating the initial liquid 
level in the tank prior to expulsion, and to enable observation of any physical processes 
occurring in the tank. Lighting of the tank interior was accomplished by bulbs mounted 
on the inner surface of the tanx wall. 


Pressurant Gas Injector Geometry 

A hemisphere injector (fig. 5) was used for all tests reported herein. The hemi- 
sphere injector was selected because it injects the pressurant uniformly in all directions 
into the ullage volume (minimizes ullage gas mixing). The uniform diffusion of pressur- 
ant is a basic assumption of the analysis used in this report. The open exit area for the 

o 

injector was 176. 8 square centimeters (27. 4 in. ). 
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Figure 4. - View port and television camera installation on test tank. Wall 
thickness, 0.762 centimeter (0.30 in.). 


Denotes pressurant gas flow paths 


Gaseous helium 



Figure 5. - Injector geometry for hemisphere injector. 
Open area, 176.8 square centimeters (27.4 in. 2 ). 

All dimensions are in centimeters (in. ). 
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Internal Tank Instrumentation 


Temperatures . - Ullage gas temperatures together with gas concentration measure- 
ments were used to determine the mass and energy content of the tank ullage. The tem- 
peratures must be obtained with sensors capable of accurate measurement of rapid 
changes in temperature. Internal tank instrumentation is illustrated in figure 6. Loca- 
tion of the vertical and horizontal ullage gas temperature rakes is indicated. The ther- 
mopile was the basic temperature measurement technique used in this investigation. The 
use of thermopiles to measure ullage gas temperature was first developed in reference 6, 
and the technique was used with good results in references 3 to 5. The main advantage of 
using thermopiles is their fast response time (between 0. 2 and 1. 0 sec) in going from 
saturated liquid to vapor. This time response is approximately an order of magnitude 
less than carbon or platinum resistance sensors which are generally used in this type of 
investigation. 

A typical three-element thermopile unit and its associated wiring schematic are il- 
lustrated in figure 7. The thermopile units were constructed of 0.202 millimeter (0.008 
in. ) chromel-constantan wire. Vertical ullage gas temperature profiles were obtained by 
stacking the individual thermopile units as shown in figure 7(b). The support structure 
was made of laminated thermoplastic to minimize thermal conduction from one measure- 
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Constantan 

Chromel 

Copper 


|Typical 



(a) Wiring schematic. 




(b) Sketch of rake. 


2.54 cm (1.0 in.) 


2.54 cm (1.0 in.) (typical 
except for lower 10 units, 
which were 1. 27 cm (0.5 in.)) 


CD -8382 


Figure 7. - Thermopile rake. 


ment station to another. The spacing between the reference and measuring levels was 
2.54 centimeters (1.0 in. ) for the top 45 thermopiles composing the vertical rake. The 
10 units at the bottom of the rake had spacings of 1 0 27 centimeters (0. 50 in. ) to obtain a 
more accurate temperature profile of the ullage gas near the liquid surface at the end of 
an expulsion. 

Platinum resistance sensors, which were located at least every tenth station starting 
from the bottom of the rake, sensed the absolute temperature at their location and pro- 
vided a reference for the thermopiles above the location. 

The horizontal rakes were composed of platinum resistance sensors spaced a maxi- 
mum of 12. 70 centimeters (5. 00 in. ) apart in the radial direction. Two platinum resist- 
ance sensors were used at each location to measure liquid and/or gas temperatures for 
the ranges of 20 to 38. 9 K (36° to 70° R) and 38. 9 to 277. 8 K (70° to 500° R). These dual 
sensors permitted more accurate measurement of liquid and gas temperatures than could 
be achieved with one sensor covering the entire range. 

The initial static temperature profile near the liquid surface was determined by a 
fixed interface rake. This rake contained 13 platinum resistance sensors spaced 0. 64 
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centimeter (0.25 in. ) apart. The range of these sensors was 20 to 38. 9 K (36° to 70° R). 
For complete tank expulsions (i.e. , 5 percent ullage to 5 percent outage), the initial 
liquid level was always assured to be on the interface rake in the thick walled tank by 
monitoring a closed- circuit television system which viewed, through a viewport in the 
tank, the position of the liquid- gas interface on a scale fixed to the interface rake (fig. 6). 
For tests conducted at the 28-, 55-, and 75-percent ullage values, the scale was moved 
lower in the tank and the television camera was refocused. Liquid level for the thin 
walled tank was determined by hot-wire liquid level sensors located on the interface rake. 

Platinum resistance sensors were also used to determine tank wall temperatures at 
12 locations and the liquid temperature at the flowmeter. Copper- constantan thermo- 
couples were used to determine tank lid temperatures at 5 locations and the pressurant 
gas inlet temperature. 

Concentrations . - The concentration of helium and hydrogen gas at five positions in 
the tank (fig. 6) was obtained by a gas sampling and analyzer system. A general sche- 
matic of this system is shown in figure 8. The sampling tubes had a 0. 157-centimeter 
(0. 062-in. ) outside diameter with a wall thickness of 0. 030 centimeter (0. 012 in. ). The 
analyzing sequencing of the five sampling tubes was done automatically. The sampling 
operation for one tube was as follows: The sample tube was purged with a negligible flow 
of helium until the liquid surface passed the entrance of the tube at which time it was in- 
cluded in the sampling sequence of the tubes above it. The gas sample was vented (also 
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at a negligible rate) to the atmosphere until its particular level was to be analyzed. At 
this time the flow control vent valve was closed and the sample flow was directed through 
a variable-leak flow control valve to a constant pressure manifold from where a portion 
of the sample passed through a 2245 cubic centimeter per hour standard leak to a mass 
spectrometer. The mass spectrometer then measured the partial pressure of the hydro- 
gen of the gas sample which was used to determine the hydrogen concentration at that 
level. The manifold and tube were then evacuated and the next level in the sequence was 
sampled and analyzed. The total time required to sample and analyze each level was 
approximately 8 seconds. 

All measurements except the gas concentrations were recorded on a high-speed dig- 

3 

ital data system. The measurements were recorded at a rate of 3. 125x10 channels per 
second. Each measurement channel was sampled every 0. 064 second. The gas concen- 
tration measurements were continuously recorded on a direct reading oscillograph. 


PROCEDURE 

The spherical test tank was filled from the bottom to approximately a 2 -percent ul- 
lage condition. It was then topped off as necessary while the tank lid and peripheral sup- 
port hardware reached steady- state operating temperatures. 

Temperature conditioning of the pressurant gas was then started. Gas flow was es- 
tablished through the heat exchanger loop, through the control valves and orifice arrange- 
ment, and then into the tank ullage from where it was vented through the conditioning line 
to the outside as shown in figure 2. The temperature control circuit shown in figure 2 
was used to get the desired pressurant gas temperature level during the flow period. 

When the pressurant gas temperature conditioning was almost completed, the liquid level 
in the test tank was adjusted to a desired value by either topping or slow draining. 

The pressurant gas flow was then stopped and the test tank was vented in preparation 
for an expulsion run. The automatic controllers and timers were preset with all the de- 
sired run and operating conditions (i. e. , tank pressure level, length of ramp period, 
length of hold period, liquid outflow valve position, start and stop times of the data re- 
cording equipment, etc. ). 

After starting the data recording equipment, the next step of the completely auto- 
matic run sequence took electrical calibrations on all pressure transducers. Immedi- 
ately following this, the test tank was pressurized over a predetermined time period to 
the nominal operating pressure of 34. 47xl0 4 newtons per square meter (50 psia). Tank 
pressure was held constant for about 30 seconds to stabilize internal temperatures. The 
tank expulsion period was then started. Approximately 90 percent of the total tank vol- 
ume was expelled at a constant volumetric flow rate. The expulsion period was stopped 
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when a hot-wire liquid level sensor located at the 95-percent ullage level indicated the 
presence of gas. The last step of the automatic rim sequence was the stopping of all data 
recording equipment. The test tank was then vented and refilled with liquid hydrogen for 
the next expulsion. 

Additional ramp pressurization runs, with no expulsion, were made for three differ- 
ent tank ullage levels. 


DATA REDUCTION 

The data reduction procedure used in this report is essentially the same that was 
used in reference 4 with the exception of the treatment for the two component ullage. The 
complete data reduction procedure is presented in appendix D. The relations derived in 
appendix D that are used in the presentation of experimental results are now summarized. 


Pressurant Gas Added (Mq) 


The actual pressurant . added to the tank ullage was determined from the relation 


M 


G, i— f 




( 1 ) 


Ideal Pressurant Requirement (Mj) 


This report uses two different relations to define the ideal pressurant requirement. 
One is used to determine the ideal requirement for the initial pressurization of the tank 
and is given by 


MP o V o 

zrt g 



( 2 ) 


This expression assumes a one component ullage. It also assumes that the pressur- 
izing gas does not mix or exchange heat with the gas in the ullage volume. For the tests 
conducted herein, the initial ullage volume contained a mixture of hydrogen and helium. 
However, because of the method used to precondition the inlet gas temperature (described 
in the PROCEDURE section), the initial ullage volume was predominantly helium. 
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Therefore, the ideal pressurant requirement used in this report is based on a 100- 
percent helium ullage (i. e. , M = 4. 003, y - 1. 67), whereas the ideal requirement used in 
reference 4 is based on 100- percent hydrogen (i. e. , M = 2. 016, y = 1. 40). In any case, 
the ideal requirement is just used as a normalizing factor and does not necessarily indi- 
cate a true minimum requirement. 

The other relation, which is used to determine the ideal requirement during the ex- 
pulsion period, is given by 


MP AV tt 

Mr - - 

1 zrt g 


(3) 


The assumption used to obtain this relation is that the incoming pressurant does not 
exchange heat or mass with the surroundings. 


Energy Balance 

Applying the first law of thermodynamics to the entire tank system (tank + ullage gas 
+ liquid) results in the expression 


£ l (“ G h G + «> dt = £ s <“L h L * + dU L> + /* f dU u + £ f 


dU. 


w 


✓ 

Total energy 
added 
(AU T ) 








Total change in liquid Total change Total change 

in tank + liquid in ullage in wall 

expelled energy energy 


(AU t ) 


(AUy) 


energy 

(AU ) 
v w' 


(4) 


Dividing equation (4) by AUrp gives 

1 _ AU L , AU n ^ AU w 
AU t AU t AU t 


(5) 


These ratios show the relative distribution of the total energy input into the system. 
The data for energy distribution presented herein are in the form of these ratios. 
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TABLE I. - MASS BALANCE RES13 

| Tank pressure, ' 

Tank Ramp 

cycle 

time Initial Mass Mass 

gec ullage added transfer 

mass, during during 

kg ramp, ramp, 

kg kg 


83 

0 .' 

762 

263 

473 

15. 

6 

22.2 

391. 

0 

428. 

8 

10.093 

0 

063 

0 . 

038 

84 



266 

479 

15. 

7 

21.6 

267. 

3 

304. 

6 

.101 


072 


,036 

85 



271 

488 

16. 

0 

21.6 

136. 

8 

174. 

4 

.108 


.073 


,044 

88 



184 

331 

14. 

6 

22.1 

396. 

2 

432. 

8 

.082 


,073 


,027 

92 



170 

306 

14. 

.5 

22.3 

138. 

0 

174. 

8 

.108 


,068 


,036 

94 



324 

583 

15. 

,9 

20.5 

395. 

4 

431. 

8 

.090 


.054 


.028 

95 



311 

560 

15. 

,7 

20.8 

271. 

2 

307. 

6 

.087 


.050 


.020 

97 


r 

335 

603 

16. 

6 

20.0 

133. 

9 

170. 

5 

.101 


.064 


.023 

7 

0 . 

409 

287 

517 

25. 

,2 

27.4 

400. 

9 

453. 

4 

0.091 

0 

.079 

0, 

.059 

10 



292 

526 

25, 

. 1 

27.6 

140. 

.5 

193. 

,2 

.086 


.073 


.051 

12 



285 

513 

24. 

.9 

27.4 

340. 

9 

393. 

2 

.085 


.071 


.047 

14 



288 

518 

24, 

.5 

28.4 

284. 

,2 

337. 

,0 

.084 


.068 


.029 

15 

1 


291 

524 

24, 

.8 

27.7 

197. 

,2 

249. 

,7 

.087 


.068 


.041 


a ±Value is probable error associated with each measurement. 


Run 

Tank 

Inlet gas 

Ramp 

Hold 

Expul- 


wall 

temper- 

time, 

time, 

Sion 


thick- 

ness, 

cm 

ature 
K °R 

sec 

sec 

time, 

sec 



.114 .034 -.013 



TABLE II. - ENERGY BALANCE RESULTS 


Run 

Ta 

wa 

thic 

nes 

nk 

ill 

;k- 

5S, 

n 

Inlet gas 
temper- 
ature 

Expul- 

sion 

time, 

sec 




Expulsion 





Energy added by 
pressurant gas, J 

Energy 
added by 
environ- 
ment, 

J 

Energy gained by tank 
wail > Au w > j 

Energy 
gained by 
ullage, 

AU U’ 

J 

Energy gained by liquid, 

AU l , j 

K 

°R 


CI 


Experimental 

(a) 

Predicted 

Experimental 

Predicted 

Experimental 

(a) 

Predicted 

83 

0.762 

263 

473 

391.0 

442.6±0.2xl0 4 

488.0X10 4 

28.6X10 4 

256.5X10 4 

259. OxlO 4 

105. 2xl0 4 

94. 4*43. 3xl0 4 

98.3X10 4 

84 



266 

479 

267.3 

414. 8±0. 2 

431.0 

19.6 

228.3 

232.8 

118.8 

93. 7 ±36 . 8 

67.8 

85 



271 

488 

136.8 

350. 2±0. 2 

355.7 

9.9 

179.9 

175.7 

112.5 

46. 7 ±21. 4 

34.7 

88 



184 

331 

396.2 

345. 5±0. 2 

401.9 

29.0 

175.1 

174.8 

157.0 

77. 0±44. 5 

99.7 

92 



170 

306 

138.0 

278. 1±0.2 

296.3 

10.1 

108.7 

106.4 

125.0 

86. 6±19. 9 

34.8 

94 



324 

583 

395.4 

529.0±0.2 

593.6 

29.0 

313.5 

354.9 

117.3 

99. 2±39. 0 

130.0 

95 



311 

560 

271.2 

502. 1±0. 3 

511.1 

19.7 

292.0 

312.7 

117.2 

66. 8±34. 6 

68.2 

97 

■ 

I 

335 

603 

133.9 

411. 3±0. 2 

401.6 

9.8 

204.7 

227.8 

100.3 

142. 0±14. 3 

33.8 

7 



517 


433. 1±0. 2xl0 4 

466.9X10 4 

30.0X10 4 

241.2X10 4 

254.9X10 4 

115. lxlO 4 

85. 6*28. 8xl0 4 

105.6X10 4 


■ 

■ 


526 

140.5 

355.5*0.1 

356.6 

10.3 

171.2 

157.0 


89.7*28.6 

36.7 

12 


■ 


513 


428. 7*0. 1 

437.9 

24.9 

233.8 

220.9 

113.1 


89.1 

14 




518 

284.2 

417. 1*0. 1 

423.3 

20.8 

223. 1 


114.5 

119. 1*27.5 

74.9 

15 

■ 

■ 

291 

524 

197.2 

381.1*0.1 


14.4 

195. 1 

188.2 

112.3 

86. 9±22. 9 

51.8 


a ± Values is probable error associated with each measurement (absolute value). 










Error Analysis 


An analysis was performed to determine the magnitude of probable error which could 
be present in the integrations of equations (D4), (D13), and (D14). Probable error is de- 
fined as follows: There is a 50-percent probability that the error will be no larger than 
the value stated. This analysis considered the errors introduced by the inaccuracies of 
temperature transducers as well as the tank pressure sensor. These calculations were 
performed for all runs for the expulsion period. The results of this analysis are included 
with the tabular expulsion data in tables I and II. No error analysis was performed on 
parameters which were dependent on the measured gas concentration data. The actual 
uncertainty of determining the gas concentration was unknown although it was estimated 
to be between ±10 to ±30 percent. 


RESULTS AND DISCUSSION 

The main parameter used to compare the effectiveness of the various operating pa- 
rameters on the amount of pressurant gas used is the nondimens ional ratio Mj/Mq. The 
symbol Mj is defined as the ideal helium pressurant mass required to pressurize or ex- 
pel a given volume of liquid at a given inlet gas temperature and tank pressure assuming 
no heat or mass transfer, and Mq is the actual pressurant requirement for the same 
inlet gas temperature and tank pressure „ A high Mj/Mq ratio means less energy and 
mass exchange. It does not necessarily mean a lower absolute pressurant requirement 
Mq as will be illustrated later in this section. 

A value of Mj/Mq equal to 1 implies that there is no heat transferred to the tank 
wall or liquid and no mass transfer. This means that for no environmental heating the 
terms AU w /AU,p and AUp/AUrp in equation (5) of the DATA REDUCTION section are zero 
and that AU^/AU T is equal to 1; that is, all the energy (AU T ) added to the tank during 
expulsion appears as an increase in ullage energy (AUjj). Therefore, any value of 
Mj/Mq or AU u /AU t less than 1 means energy is lost by the ullage system. This loss 
of ullage energy would then appear as a change in tank wall energy and/or liquid energy; 
that is, AU w /AU t and/or AUp/AU,p would be greater than zero. 

The discussion of results first presents the effects of the various operating param- 
eters on the ratio Mj/Mq for the expulsion period only, which was of primary interest, 
followed by the mass transfer Mj/Mq results. Then the results of the energy balances 
are presented in an attempt to point out major reasons for the ratios Mj/Mq or 
AUjj/AU^ being less than 1. Finally, a comparison is made between the experimental 
results and the analytically predicted results to determine the validity of the analytical 
program. The analytical results are presented in the figures along with the correspond- 
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Inlet gas 
temperature, 
K (°R) 


Runs 


□ 322(580) 94, 95, 97 

A 267 ( 480.6) 83, 84, 85 

O 183 (329) 88, 92 



2l I I I I I I 

100 150 200 250 300 350 400 


Total expulsion time, sec 


Figure 9. - Comparison of ideal pressurant requirement to actual 
pressurant requirement ratio as function of expulsion time for 
three inlet gas temperatures. Tank pressure, 34.47x104 newtons 
per square meter (50 psia); wall thickness, 0. 762 centimeter 
(0.30 in.). 


ing experimental results. The comparison between experimental and analytical results 
are given in terms of an average deviation, which is defined as 


N 



(Experimental rati o) - (Analytical ratio) 
(Experimental ratio) 




( 100 ) 


(27) 


where N is the number of data points in a given set of operating conditions (i. e. , for a 
constant inlet gas temperature of 322 K (580° R), N would be 3 for the data presented in 
fig. 9). For convenience, all deviations between the experimental and analytical results 
are summarized in table in. 

The operating parameters (e. g. , inlet gas temperature, outflow rate, and tank wall 
thickness) and major experimental and analytical results are summarized in tables I 
and II. Table I gives the experimental and analytical mass balance results, while table II 
gives the corresponding energy balance results. 


Effect of Inlet Gas Temperature 

Pressurant requirements . - The effect of inlet gas temperature is shown in figure 9 
on the basis of Mj/Mq for various expulsion times. Expulsion time is the total time 
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TABLE m. - DEVIATIONS BETWEEN EXPERIMENTAL AND ANALYTICAL 


RESULTS FOR HEMISPHERE INJECTORS 


Run 

Inlet gas temperature 

Tank wall 

Percent deviation between experimental 


K 

°R 

thickness, 

and analytical results a 


cm 

Mj/Mq 

au w /au t 

AU l /AU t 

83 

263 

473 

0.762 

+10.25 

-0.91 

-4.00 

84 

266 

479 



+3.92 

-2.09 

+27.77 

85 

271 

488 



+ 1.57 

+2.20 

+26. 15 

88 

184 

331 



+16.38 

0 

-29. 12 

92 

170 

306 



+6.54 

+2.11 

+59.80 

94 

324 

583 



+12.21 

-13.16 

-30.89 

95 

311 

560 



+ 1.79 

-7. 15 

-2.34 

97 

335 

603 



-2.40 

-11.31 

+76.26 

7 

287 

517 

0.409 

+ 7.82 

-5.75 

-23.24 

10 

292 

526 



+ .25 

+8.33 

+59. 18 

12 

285 

513 



+2.14 

+5.43 

+4.39 

14 

288 

518 



+ 1.48 

+5.88 

+37. 13 

15 

291 

524 



+ .84 

+3.60 

+40.36 


a Analysis overpredicts, analysis underpredicts, +. 


required to expel liquid from a 5-percent ullage to a 95-percent ullage. Therefore, each 
data point represents a complete expulsion. For a given inlet gas temperature, there is 
an increasing pressurant requirement (decreasing Mj/Mq) for increasing expulsion 
times. The longer the pressurant (ullage) gas is exposed to the cold surroundings, the 
greater the loss in pressurant energy. The ratio Mj/Mq also decreases for increasing 
inlet gas temperature (from 0. 540 for the 183 K (329° R) inlet gas temperature at a 138- 
second expulsion to 0. 263 for the 322 K (580° R) inlet gas temperature at a 395-second 
expulsion). This decreasing ratio implies that a larger percentage of the pressurant gas 
energy is lost to the tank wall and/ or to the liquid as the inlet gas temperature is increased. 
The values of Mj/Mq obtained here using helium as the pressurant are approximately 
the same as those obtained in reference 4 using hydrogen as the pressurant for similar 
inlet gas temperatures and expulsion times. A comparison of the actual pressurant re- 
quirements M g for the three inlet gas temperatures for various expulsion times is 
shown in figure 10. The actual pressurant requirements (Mq) decrease for increasing 
inlet gas temperature and increase for increasing expulsion times. There is an average 
decrease of 1. 3 percent in Mq for a 10 K (18° R) increase in inlet gas temperature. 

The results of reference 4 indicated that there was an average decrease in pressurant 
gas requirements of 1. 4 percent for a 10 K (18° R) increase in inlet gas temperature 
when using hydrogen as the pressurant gas. 
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Open and closed symbols denote experimental 
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Figure 10. - Comparison of actual pressurant added as function of 
expulsion time for three inlet gas temperatures. Tank pressure, 
34.47X10 4 newtons per square meter (50 psia) ; wall thickness, 

0. 762 centimeter (0. 30 in.). 


Although the absolute value of Mq decreases for increasing inlet gas temperatures, 
the ratio Mj/Mq at 322 K (580° R) is less than the Mj/Mq ratio at 183 K (329° R) be- 
cause of the greater decrease in ideal requirements. 

The shaded symbols shown in figures 9 and 10 are the results as predicted by the 
analytical program. The agreement between the analysis and experimental results is 
better for the 322 and 267 K (580° and 480. 6° R) inlet gas temperatures (see table III for 
actual deviations) than for the 183 K (329° R) inlet temperature. The reason for the 
greater deviation for the lower temperature is not clearly understood. The deviation be- 
tween the analytical and experimental results is also greatest for the longest expulsion 
time in all cases. This could be attributed to the diffusion of hydrogen into the helium 
ullage (to satisfy its partial pressure) which is a time- dependent process. The longer 
the liquid hydrogen is exposed to a helium ullage the more hydrogen will diffuse into the 
ullage. The limiting condition is when the partial pressure of hydrogen in the ullage is 
satisfied. The addition of hydrogen into the ullage would therefore reduce the actual 
helium pressurant requirement. The analysis does not consider this mass transfer and 
would therefore overpredict the pressurant requirements, especially for the longer ex- 
pulsion times. 

Mass transfer . - The amount of mass transfer was not directly measured experi- 
mentally. It was determined indirectly by the use of equation (D8). The accuracy of de- 
termining the mass transfer is strongly influenced by the accuracy in which the gas sam- 
pling and analysis system determines the helium-hydrogen concentrations in the ullage. 

A comparison of the ratio of mass transferred during expulsion to the actual pres- 
surant added to the tank M^/Mq is presented in figure 11 for different expulsion times 
and the three inlet gas temperatures. The experimental results presented in this figure 
indicate that helium is absorbed into the liquid hydrogen for the short expulsion times for 
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Runs 


Inlet gas 
temperature, 
K(°R) 



Figure 11. - Comparison of mass transfer to mass added ratio as function of expul- 
sion time for three inlet gas temperatures. Tank pressure, 34.47x10^ newtons 
per square meter (50 psia); wall thickness, 0.762 centimeter (0.30 in.). 


the three inlet gas temperatures. Both the 322 and 267 K (580° and 480. 6° R) inlet gas 
temperatures show a trend toward evaporation of liquid for increased expulsion times. 
This trend is not apparent for the 183 K (329° R) inlet gas temperature. 

As mentioned previously, the determination of mass transfer is strongly influenced 
by the determination of the helium -hydrogen gas concentration in the ullage. Figure 12 
compares the ullage gas concentration profiles prior to and after expulsion for a run 
which indicated absorption (run 85) and a run which indicated evaporation (run 83). Both 
runs show an increase of hydrogen gas in the ullage over the expulsion period. This fact 
alone indicates evaporation and/or diffusion. However, the percent of helium near the 
liquid surface (where most of the ullage mass is concentrated) for run 85 was much less 
than that of run 83, resulting in a much lower final ullage mass. The use of equation (8) 
(appendix D) for run 85 results in a value of mass transfer (absorption) of 0. 458 kilogram 
(1. 009 lb). Based on a 1. 2-mole-percent solubility of helium in liquid hydrogen (ref. 7) 
the maximum amount of helium that could go into solution is approximately 2. 618 kilo- 
gram (5. 770 lb). Even though all absorption values were well below this maximum value, 
the actual mass-transfer values are considered inconclusive because of the rather large 
uncertainties in determining concentration gradients. 

There are no analytical comparisons for the mass transfer because the analysis 
neglects mass transfer in its development. 

Energy remaining in ullage . - The ratios of the energy increase in the ullage over 
the expulsion period to the total energy added to the system AU^/AU ^ for different ex- 
pulsion times are compared in figure 13 for the three inlet gas temperatures. For all 
runs, between 21. 0 and 43. 5 percent of the total energy that was added to the system re- 
mains in the ullage after expulsion. In general, the ratio AU^/AUrp decreases for in- 
creasing inlet gas temperatures for a constant expulsion time. Also, for a given inlet 
gas temperature, the ratio AU^/AUrp decreases with increasing expulsion time. It 
should be noted that the absolute value of AUpj does not change significantly for the 
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Open and closed symbols denote 
runs 85 and 83, respectively 
Flagged symbol indicates value 
the same for both runs 
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Figure 12. - Comparison of ullage gas concentration 
profiles prior to and after expulsion for runs 83, 
which indicated evaporation, and 85, which indicated 
absorption. Tank pressure, 34. 47xl0 4 newtons per 
square meter (50 psia) ; inlet gas temperature: for 
run 85, 271 K (488° R); for run 83, 263 K (473° R). 
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Figure 13. - Comparison of energy increase in ullage to total energy 
added ratio as function of expulsion time for three inlet gas tem- 
peratures. Tank pressure, 34.47X10 4 newtons per square meter 
(50 psia); wall thickness, 0. 762 centimeter (0. 30 in. ). 


three operating inlet gas temperatures used during testing. The mean increase in ullage 
energy for these series of runs (table II) was 117x10^ joules (1106 Btu) with a standard 
deviation of 13. 7x10^ joules (129 Btu). Any trends in the ratio AlIy/AU-p, therefore, 
depend mainly on variations in the total energy added (AU-p) to the system due to varia- 
tions in energy losses to the tank wall and liquid. 

Energy added to tank wall . - The ratios of energy gained by the tank wall to the total 
energy added to the system AU w /AUp for different expulsion times are compared in 
figure 14 for the three inlet gas temperatures. In general, between 37. 5 and 56. 2 per- 


Inlet gas Runs 

temperature, 

K (°R) 

□ 322+22 (580+40) 94, 95, 97 

V 267+4 (480.6+7) 83, 84, 35 
O 183+11 (329+20) 88, 92 

Tailed symbol indicates experimental 



Expulsion time, sec 

Figure 14. - Comparison of energy gained by wail to total energy 
added ratio as function of expulsion time for three inlet gas tem- 
peratures. Tank pressure, 34.47X10 4 newtons per square meter 
(50 psia); wall thickness, 0. 762 centimeter (0. 30 in. ). Curves 
are faired through experimental data. 
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cent of the total energy added to the system was gained by the tank wall over the range of 
conditions. The results of reference 4 showed that between 30 and 58. 8 percent of the 
total energy added to the tank was gained by the tank wall for similar test conditions when 
using hydrogen as the pressurant. A comparison of results indicate that there is gener- 
ally a 5 to 7. 5 percent increase in the percent of total energy that is lost to the tank wall 
when using helium as the pressurant. As can be seen in table II and in figure 14, both 
the absolute value of AU W and the ratio AU w /AUrp increase with increasing inlet gas 
temperature. The increase in AU W is due to the larger driving potential AT for heat 
transfer between the ullage gas and the tank wall. The total energy added to the system 
(AUrp) does not increase in the same proportion as AU W resulting in the increased ratio 

au w /au t . 

The agreement between the analysis and experimental value of AU w /AUrp are very 
good for the 183 and 267 K (329° and 480. 6° R) inlet gas temperature (within ±2.2 per- 
cent, see table III). However, the agreement for the 322 K (580° R) inlet gas is not as 
good (within ±13.2 percent). 

Energy gained by liquid . - Figure 15 is a comparison of the ratio of energy gained 
by the liquid to the total energy added to the system AU^/AU-p for different expulsion 
times for the three inlet gas temperatures. In all cases, between 12.8 and 33.7 percent 
of the total energy added to the system appears as an increase in liquid energy (liquid 
heating). The results of reference 4 indicated that between 15 and 25 percent of the total 
energy added to the system appears as an increase in liquid energy for similar test con- 
ditions using hydrogen as the pressurant. There is a relatively large amount of scatter 
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Figure 15. - Comparison of energy gained by liquid to total energy 
added ratio as function of expulsion time for three inlet gas tem- 
peratures. Tank pressure, 34.47X10 4 newtons per square meter 
(50 psia); wall thickness, 0. 762 centimeter (0. 30 in. ). Curves 
are faired through experimental data. 
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in the data shown in figure 15. This scatter could be the result of the error in experi- 
mentally determining AU T . The probable error associated with each experimental de- 

J - J .44 

termination of AU-^ is between 10 and 58 percent (i.e., 14.3x10 out of 142.0x10 J and 

44. 5X10 4 out of 77. OxlO 4 J, see table H). 

The analytical predictions of the ratio AUjVAUrp are also presented in figure 15 
(see table in for average deviations). The large discrepancy between the analysis and ex- 
perimental results could also be the result of the error in determining AU^. 

Temperature distributions . - The results that were discussed previously point out 
that between 72. 4 and 88. 6 percent of the total energy that was added to the system is 
either absorbed by the tank wall or remains in the ullage. The correlation between the 
analysis and experimental data, therefore, depends largely on the ability of the analysis 
to predict final wall and ullage gas temperature profiles. These profiles are, in turn, 
used to determine the increase in wall and ullage energy and the final ullage mass. The 
ability to predict these temperatures explains the fair agreement between experimental 
data and analysis reported in references 1 and 3 to 5. 

Figure 16 presents comparisons of experimental and analytical ullage gas and wall 
temperature profiles for the three inlet gas temperatures. Figure 16(a) presents a com- 
parison of experimental and analytical wall and ullage gas temperature profiles at the end 
of a 396-second expulsion for an inlet gas temperature of 184 K (331° R). For this run 
the analysis accurately predicted the AU w /AUrp ratio but underpredicted the Mj/M^ 
ratio by 16.38 percent. The experimental gas temperatures shown in figure 16(a) were 
obtained from the vertical rake. The horizontal rakes indicate these are average radial 
temperatures at their respective vertical position. In the absence of any mass transfer, 
the pressurant mass required for an expulsion could be determined as the difference be- 
tween the final mass in the ullage and the initial mass prior to expulsion. 

One of the inputs to the analytical program is the initial experimental temperature 
profile prior to expulsion. The analytical program uses this profile together with the 
temperature-density relation for helium to obtain the initial ullage mass. The initial 
ullage mass as determined by the analysis is, in all cases, larger than determined ex- 
perimentally because the experimental values contain hydrogen in the ullage. However, 
the difference between the analytically and experimentally determined initial ullage mass 
represents only a small fraction (less than 3 percent) of the final ullage mass. There- 
fore, the deviation between the analytical and experimental pressurant requirements 
would largely be the result of the predicted final ullage temperature profile. As seen in 
figure 16(a) the analytical gas temperatures are slightly lower in the lower 20 percent of 
the tank but begin to get higher toward the top of the tank. The overprediction of temper- 
atures in the majority of the ullage results in a 16. 38-percent underprediction in the final 
ullage mass. 

The analytically predicted wall temperatures (fig. 16(a)) are very close to those ob- 
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(c) Expulsion time, 134 seconds; inlet gas temperature, 335 K (603° R ); run 97. 

Figure 16. - Concluded. 

served experimentally in the lower 47 percent of the tank ullage but begin to get higher 
toward the top of the tank. Even though the calculated wall temperatures are higher in 
the upper portion of the tank wall, the analysis predicts the energy gained by the tank wall 
(AU ) accurately. The effect of overpredicting the temperatures in the upper portion of 
the tank wall is offset by an underprediction of the temperature of the tank lid (point zero 
in fig. 16(a)), which is approximately 22 percent of the total tank mass. 

Figure 16(b) is a comparison of experimental and calculated wall and ullage gas tem- 
perature profiles at the end of a 137-second expulsion for an inlet gas temperature of 
271 K (488° R). Here the calculated gas temperatures agree very closely with the exper- 
imental temperatures in the lower 50 percent of the ullage but begin to get lower toward 
the top of the tank. However, the predicted pressurant requirement is only 1. 57 percent 
higher (lower Mj/Mq) than that observed experimentally. 
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The calculated wall temperatures agree fairly well with the experimental tempera- 
tures except near the top of the tank where they begin to get higher. However, this over- 
prediction is again offset by an underprediction of the lid temperature resulting in an 
under prediction for wall energy gained of 2. 20 percent. 

A comparison of experimental and calculated wall and ullage gas temperatures at the 
end of a 13 4- second expulsion for an inlet gas temperature of 335 K (603° R) is shown in 
figure 16(c). Here again there is good agreement between the calculated and experimen- 
tal gas temperatures throughout the tank ullage resulting in a good prediction for the 
pressurant requirements (within 2. 40 percent). However, in this case the lower calcu- 
lated lid temperature does not offset the higher calculated wall temperatures in the upper 
portion of the tank. And the analysis overpredicts the energy gained by the tank wall by 
11.31 percent. 

The results of these comparisons indicate that the analytical program is able to pre- 
dict with reasonable accuracy the pressurant gas requirement and the major energy 
losses for the test conditions imposed herein. The results also indicate that the analyti- 
cal assumption of a one component ullage (no mass transfer) is adequate and that what 
hydrogen gas there is in the ullage has little influence on helium pressurant requirements 
or energy distributions. 


Effect of Tank Wall Thickness on Pressurant Requirements 

The results of reference 3 indicated that between 77 and 93 percent of the total en- 
ergy lost by the pressurant gas was lost to the tank wall. It was expected, therefore, 
that tank wall thickness would play an important role in determining the pressurant re- 
quirements during the expulsion period. Additional tests were performed using a 1. 52- 
meter- (5-ft-) diameter spherical tank similar to the one already discussed except the 
tank wall was chem-milled to an average thickness of 0. 409 centimeter (0. 161 in. ). Al- 
though the chem-milling reduced the wall thickness by 46. 3 percent, the total weight of 
the tank was reduced by only 29. 5 percent because of the structural requirement of keep- 
ing the same thickness for the tank neck, lid, and girth support. 

Figure 17 presents a comparison of the Mj/Mq ratio for different expulsion times 
for the two tank wall thicknesses (0. 762 and 0. 409 cm; 0. 30 and 0. 161 in. ). Each tank 
used the hemisphere injector. The average inlet gas temperature for the 0. 409- 
centimeter (0. 161-in. ) tank was 290 K (~520° R). The average inlet gas temperature for 
the 0. 762- centimeter (0.30-in.) tank was 267 K (~480° R). The 0. 409-centimeter- 
(0. 161-in. -) thick tank has an average of 6. 3 percent less pressurant requirement, 
higher Mj/M^ ratio, than the 0. 762-centimeter- (0. 30-in. -) thick tank for all expulsion 
times. Similar results were obtained in reference 4 using hydrogen. In all cases, the 
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Figure 17. -Comparison of ideal pressurant requirement to actual 
pressurant requirement ratio as function of expulsion time for 
two wall thicknesses. Tank pressure, 34. 47x10^ newtons per 
square meter (50 psia); inlet gas temperature, 277±10 K 
(499±18° R). 


analysis overpredicts the pressurant requirements (under predicts the Mj/Mq ratio). 

See table III for deviations between analytical and experimental results. 

The 6.3-percent decrease in pressurant requirements was less than expected. How- 
ever, the test with the 0. 409- centimeter- (0. 161-in. -) thick tank had an average 23 K 
(41° R) higher inlet gas temperature. Based on the results already obtained on the effect 
of inlet gas temperature on pressurant requirements (fig. 9), this 23 K (41° R) higher in- 
let temperature would normally decrease the pressurant requirement by approximately 
2. 9 percent. Thus for equal inlet gas temperatures, a nominal decrease in pressurant 
requirements of approximately 9. 0 percent would be expected for the 0. 409- centimeter - 
(0. 161-in. -) thick tank. 

Figure 18 is a comparison of AU w /AU,p for different expulsion times for the two 
tank wall thicknesses. In this figure the AU /AU t ratio is broken down into the frac- 
tional portion gained by the tank wall and the total fraction gained by the tank wall plus 
tank lid. (The tank lid represented 30 percent of the total tank weight for the thinner wall 
tank compared to 22 percent for the thicker tank. ) The fraction of the total energy added 
that is gained by the tank wall is approximately the same for both tank thicknesses (lower 
curve of fig. 18) and ranges between 0.37 and 0. 40. The values of AU w /AU,p that in- 
clude the lid indicate a 4- percent decrease in AU w /AU t for the thinner wall tank. It is 
apparent then that the reduction in the total AU /AU™ is almost entirely due to a reduc- 
tion in the fraction of the total energy absorbed by the lid for the thinner wall tank. The 
reason for the thinner tank’s reduced lid energy gain was that its initial lid temperature, 
prior to expulsion, was an average of 41 K (73.8° R) higher than the thicker tank's. This 
higher initial lid temperature reduced the temperature difference between the lid and the 
ullage gas resulting in less energy gained. The absolute value of the total wall energy 
gained (table II) by the 0. 409-centimeter- (0. 161-in. -) thick tank was approximately 
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Figure 18. - Comparison of energy gained by wall to total energy added 
ratio for two tank wall thicknesses. Tank pressure, 34. 47x1$ 
newtons per square meter (50 psia) ; inlet gas temperature, 277+ 10 K 
(499°±18° R). 


6. 9 percent less than for the 0. 762 -centimeter- (0. 30- in. -) thick tank. 

The absolute values of increase in ullage energy, liquid heating, and mass transfer 
are approximately the same (within measurement uncertainty) for the two tank wall thick- 
nesses (see tables I and II). Thus, the 6.3-percent reduction in pressurant requirements 
when using the 0. 409- centimeter- (0. 161-in. -) thick wall tank is almost entirely due to 
the 6. 9-percent reduction in tank wall heating (AU W ). 


Pressurant Requirements for Initial Pressurization 

The amount of pressurant gas needed to initially pressurize a propellant tank may be 
important for certain missions, particularly for multiburn missions where the tank is 
vented after each burn or where the coast period between firings is long enough to enable 
the ullage gas to collapse. 

As stated in the INTRODUCTION, the purpose of this investigation was to determine 
the capability of the analysis to predict the pressurant requirements during the initial 
pressurization period as well as the expulsion period. For this purpose, data were col- 
lected during the initial pressurization period for various pressurizing rates, inlet gas 
temperatures, and ullage volumes. 

Figure 19 is a comparison of the Mj/Mq ratio as a function of inlet gas temperature 
for two ramp rates at an initial ullage volume of approximately 4 percent. The data were 
taken using the hemisphere injector. At constant inlet gas temperature, the pressurant 
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Ramp rate, 
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Figure 19. - Comparison of ideal pressurant requirement to actual pres- 
surant requirement ratio for ramp period as function of inlet gas temper- 
ature for two ramp rates. Ullage volume, 4 percent. 

requirements decreased for increased ramp rates. The Mj/Mq ratio decreases for in- 
creasing inlet gas temperature. The reduced pressurant requirements (for the ramp 
period) for fast ramp rates using low inlet gas temperature is in agreement with the 
trends during the expulsion period. 

The modification of the analysis of reference 1 for the ramp period is discussed in 
appendix C. As can be seen in figure 19, the analysis is not capable of predicting the 
pressurant requirements accurately during the initial pressurization of the 4-percent 
ullage. However, the prediction of the total pressurant requirement (initial pressuriza- 
tion and expulsion) is still adequate because the amount of gas required to pressurize the 
4-percent ullage initially was only 4. 0 percent (maximum) of the pressurant requirement 
during expulsion. The absolute pressurant requirements (both experimental and analyti- 
cal) for the data presented in figures 19 and 20 are given in table IV. 

The comparison of Mj/Mq as a function of initial ullage volume for various ramp 
rates at an inlet gas temperature of approximately 180 K (324° R) is shown in figure 20. 
This figure indicates increased Mj/Mq for increased ullage volume for any given ramp 
rate. There is an increased Mj/Mq for increased ramp rates for a given ullage volume. 
As can be seen in figure 20, the accuracy of the analytical prediction is not consistent. 

The actual deviation between analytical and experimental values range between 0 and 
20. 5 percent. The transient process that occurs during the initial pressurization of the 
tank is too complex to be described by the present analytical model. The analytical pro- 
gram can, however, at least be used to predict the approximate magnitude of pressurant 
requirements during the ramp period. 
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Figure 20. - Comparison of ideal pressurant requirement to actual pres- 
surant requirement ratio for ramp period as function of initial ullage 
volume for three ramp rates, inlet gas temperature, 180 K (324° R). 
Curves are faired through experimental data. 


TABLE IV. - COMPARISON OF EXPERIMENTAL AND ANALYTICAL VALUES 


OF PRESSURANT GAS REQUIREMENTS FOR RAMP PERIOD 


Run 

Initial 

Ramp rate 


Mass added 


Inlet gas 


ullage 


o 





temperature 


volume, 

psi/sec 

(N/na )/sec 

Experimental 

Analytical 




percent 



lb 

kg 

lb 

kg 

K 

°R 

84 


1 

1.65 

11.37X10 3 

0. 159 

0.072 

0. 148 

0.067 

265 

478 

86 



.73 

5.03 

. 177 

.080 

.202 

.092 

209 

377 

88 



1.65 

11.37 

. 160 

.073 

.247 

. 112 

178 

320 

95 



1.65 

11.37 

. 110 

.050 

.175 

.079 

294 

530 

96 



. 73 

5.03 

.101 

.046 

. 173 

.078 

302 

545 

91 

\ 


.86 

5.92 

.197 

.089 

.247 

.112 

177 

319 

124 

28 

.86 

5.92 

1.344 

.610 

1.179 

.535 

181 

325 

125 

28 

.47 

3.24 

1.811 

.821 

1.701 

.772 

181 

326 

147 

55 

.86 

5.92 

2.245 

1.018 

2.312 

1.049 

174 

313 

148 

55 

.73 

5.03 

2.360 

1.070 

2.389 

1.084 

180 

323 

149 

55 

.47 

3.24 

2.257 

1.024 

2.674 

1.213 

182 

328 

394 | 

75 

.86 

5.92 

2.750 

1.248 

3.054 

1.385 

169 

304 

395 

75 

. 73 

5.03 

2.618 

1. 188 

2.350 

1.066 

166 

299 

396 

75 

.47 

3.24 

2.522 

1. 144 

2.535 

1. 150 

168 

302 


35 



SUMMARY OF RESULTS 


Tank pressurization and propellant expulsion tests were conducted in 1.52-meter- 
(5 -ft-) diameter spherical tanks to (1) determine what factors have the greatest influence 
on pressurant gas requirements when using helium as the pressurant, and (2) verify the 
capability of an analysis to predict the helium pressurant requirements during the initial 
pressurization period as well as the expulsion period. Tests were conducted using three 
inlet gas temperatures over a range of liquid outflow rates. The results of this investi- 
gation are now given. 


Experimental Results 

The experimental results indicate an average decrease in helium pressurant gas re- 
quirements (M q ) of 1.3 percent for a 10 K (18° R) increase in inlet gas temperature. 

Increased inlet gas temperature decreases the residual mass and slightly decreases 
the energy remaining in the ullage volume after a given expulsion. Of the total energy 
added to the tank system, between 37. 5 and 56. 2 percent was lost to the tank wall (of this 
between 13 and 20 percent was lost to the tank lid) and between 12. 8 and 33. 7 percent was 
lost to the liquid. 

The mass transfer values obtained in this investigation are inconclusive because of 
the uncertainty in measuring the ullage gas concentration gradients. Additional work is 
necessary to develop better techniques for measuring gas concentration gradients and 
mass transfer. 

Decreasing the tank mass by 29. 5 percent (by decreasing the tank wall thickness 
from 0. 762 to 0. 409 cm; 0. 30 to 0. 161 in. ) decreased the pressurant requirement by an 
average of 6. 3 percent (9. 0 percent when corrected for variations in in'et gas tempera- 
ture). 

The effects of inlet gas temperature and ramp rate on the pressurant gas required 
during the initial pressurization of the tank were as follows: 

1. Increased inlet gas temperature decreased the pressurant requirement Mq and 
mass ratio Mj/Mq for constant ramp rates. 

2. Increasing ramp rate decreased the pressurant requirement and increased the 
mass ratio for constant inlet gas temperatures. 

3. Increased pressurant requirement and increased mass ratio were obtained for 
larger initial ullage volumes for a constant ramp rate. 

The trends shown herein for various inlet gas temperatures, liquid outflow rates, 
and tank wall thickness are consistent with the results obtained using hydrogen as the 
pressurant (ref. 4). A comparison of the results obtained for the two pressurants indi- 
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cate that the mass ratios Mj/Mq for similar test conditions are approximately the same. 
It therefore would require approximately twice as much mass to expel liquid hydrogen 
using helium as the pressurant instead of hydrogen. 


Comparison of Analytical and Experimental Results 

The comparison between the analytical and experimental results indicate that, for 
the range of test conditions used, the analytical program and assumptions are adequate to 
allow prediction of the pressurant gas requirements during the initial pressurization as 
well as the expulsion period when using helium as the pressurant. The general results 
as predicted by the analysis are as follows: 

1. The pressurant requirements were predicted to within an average of 4. 41 percent 
for the 322 and 267 K (580° and 480.6° R) inlet gas temperatures and within 11. 46 per- 
cent for the 183 K (329° R) inlet gas temperature. 

2. Tank wall heating was predicted to within an average of 4.78 percent for all runs. 

3. Liquid heating was predicted to within an average of 33. 28 percent for all runs. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 22, 1969, 

180-31. 
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APPENDIX A 


INCORPORATION OF A VARIABLE GEOMETRY TO 
THE ANALYSIS OF REFERENCE 1 

The basic analysis used in this report to predict pressurant gas requirements was 
developed by W. H. Roudebush in reference 1 for a cylindrical tank. 

The major assumptions in the analysis of reference 1 are as follows: 

(1) The ullage gas is nonviscous. 

(2) The ullage gas velocity is parallel to the tank axis and does not vary radially or 
circumferentially. 

(3) The tank pressure does not vary spatially. 

(4) The ullage gas temperature does not vary radially or circumferentially. 

(5) The tank wall temperature does not vary radially or circumferentially. 

(6) There is no axial heat conduction in either the gas or the wall. 

(7) There is no mass transfer (condensation or evaporation). 

(8) There is no heat transfer from the pressurant gas to the liquid. 

Experiments performed at Lewis (ref. 3) confirmed most of these assumptions when 

using a diffuser type injector such as the one used herein. The experimental results in- 
dicated, however, that there is significant heat transfer from the gas to the liquid with 
resulting mass transfer. 

For the purposes of this report, the analysis of reference 1 was modified for appli- 
cation to arbitrary symmetric tank shapes, and an attempt was made to incorporate the 
heat transfer from the gas to the liquid. The treatment of internal hardware (e.g. , tank 
baffles, instrumentation) was also modified to correspond to the treatment of heat trans- 
fer to the tank wall. 

The primary equations which deal with the pressurizing gas upon entering the tank 

are: 

(1) The energy equation 

(2) The continuity equation 

(3) The equation of heat transfer for a point in the tank wall 


v 

ENERGY EQUATION 

The form of the energy equation used in the analysis in reference 1 for cylindrical 
tanks is 
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Modifying this equation to account for both arbitrary symmetric tank shapes and internal 
tank heat sinks gives 
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The first term on the right includes the effect of wall curvature. The last term, the 
energy lost to the internal hardware, is treated as the summation of hardware compo- 
nents: (1) laminated thermoplastic, (2) stainless steel, and (3) copper. For the tanks in 
this investigation, 

% _ a hM t h " T G> , a9 

L p H v„ ‘ 

Hardware 

components 


Gluck and Kline, in reference 8, employed the free convection correlation to the 
pressurant gas (hydrogen, helium) for the pressurized transfer of liquid hydrogen: 


= Nu = 0. 13 (GrPr) 1 / 3 (A3) 

k 


This correlation is used herein even though it was developed for cylindrical tanks. 
Pressurant gas transport properties were evaluated at the mean of the gas and wall tem- 
peratures. 


CONTINUITY EQUATION (AREA = f(x)) 

The basic form of the continuity equation for a cylindrical tank is presented in refer- 
ence 1 (eq. (24)) as 


dV = ^l_ / 3T + V 9T\ Z 2 3P 
ax zt Vat dx ) zp at 
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The modified form of the continuity equation due to variations in tank radius with 
distance along the vertical axis becomes 

3V _ Z 1 / 0T + V 3T \ Z 2 3P 2V 3r 
3x ZT \ 3t 3 x / Zp 3t r dx 


where Z^ and Zg are defined in reference 1 as 



The last term in equation (A4) evolves from the derivation as follows. For the one- 
dimensional expression for continuity, 

— (pV A) + — (pA) = 0 
0X 0t 

2 

The substitution A - nr is made where r is the position radius at location x along 
the vertical axis: 


— (pV r 2 ) + — (pr 2 ) = 0 
dx 3t 


The expression for density from the equation of state p = MP/ZRT is substituted: 




= 0 


The following velocity equation is obtained after performing the partial differentiation and 
after rearranging terms: 


3V = 

1 + 1 f dZ ) 

fdT V dT) 

1 f dZ ) - 1 

3x 

T Z \3T/p 

\ dt dx J 

Z \3P Jrp P 


When the expressions involving Z^ and Z 2 are substituted in 
tion (A4) is obtained. 


3P _ 2 V dr 

3t r dx 

this equation, equa- 
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TANK WALL HEAT TRANSFER 


Reference 1 (eq. (18)) gives the heat-transfer equation which represents the change 
in wall temperature as a result of the convective process for a cylindrical tank: 

— = ^ (T - T ) + (A5) 

3t *w*>w C w Vw C w 

where q w is the rate of heat addition per unit area to the tank wall from outside the tank. 
For a small element of volume in the x-direction, equation (A5) can be written as 

St 

P "' C » V lf * h c A(T -V + «w < A6 > 

For a wall of arbitrary shape, the following is evident from the sketch: 



A _ 2jtt As _ As _ _1_ 

V 27rr Ar Ax Ar Ax l 

Therefore, equation (A5) holds also for this case. 

To account for the large mass concentration at the top of the tank, an equivalent l w 
was used. This l was obtained by dividing the mass of the tank lid and flange connec- 
tions by the surface area at the first net point. 
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APPENDIX B 


EQUATIONS OF HEAT AND MASS TRANSFER AT THE GAS-LIQUID INTERFACE 

The energy and continuity equations (Al) and (A4) should be modified to incorporate 
both heat transfer from the ullage gas to the liquid surface and mass transfer into the 
analysis. 

The energy equation should incorporate three additional terms: 

(1) The heat-transfer rate from the ullage gas to the liquid interface (q-^^g) 

(2) The heat- transfer rate from the interface to the liquid (qg_ L ) 

(3) The energy associated with mass transfer (M^X) (evaporation is positive) 

These additional terms can be related by performing an energy balance at the gas- liquid 
interface as done by W. A. Olsen in reference 9. The resulting energy relation is based 
on the assumption requiring the interface for a pure system to be at the saturation tem- 
perature corresponding to the tank pressure. As shown by the sketch, 


Interface 



the energy balance at the interface is given by 

M t X 

q U-S - q S— L + 

In the case where helium is the pressurant gas for liquid hydrogen, it is assumed 
that there is sufficient evaporation to maintain a hydrogen blanket over the entire liquid 
surface. The interface can then be assumed to be at the saturation temperature for a 
pure hydrogen system. 

Experimental data indicate die energy associated with the mass transfer was rela- 
tively small so that the following assumption is made: 



With the further qualification that the environmental transfer to the liquid is also small, 
it follows that 


Vs a Vl “ £ ®L> 


(B2) 


J 

The term — (U^) can be determined from experimental data. The analysis requires 
. ^ 

that q^g must be related to the ullage gas variables. This is done by the following 
equation which involves defining a heat-transfer coefficient (h c j ) and a temperature 
(Tg) somewhere in the vapor: 

Vs * »c, l< t S - T sat> < B3 > 


The flow process is free convection flow of pressurant gas essentially down the tank wall 
and then radially inward across the liquid surface. 

With regard to h c reference 10 develops an equation from boundary layer theory 
for forced flow across a horizontal, semi-infinite, constant temperature flat plate given 

by 


Nu 






(B4) 


The velocity of the gas across the liquid surface in terms of the gas velocity Vq 
down a vertical wall is given in reference 11 as 

V L = 0. 0975 V G (B5) 

where Vq, obtained by solving the integrated energy and momentum equations at the wall 
boundary, is given by 


= 1. 185 — 

0 pz 


Gr 


1/2 


1 1/2 


(B6) 


(l+ 0.494 Pr 2 / 3 ) 1 

Combining equations (B5) and (B6) for and substituting into (B4) gives 

Nu = Vi = 0-226 Pr 1 / . 3 G r 1 / 4 = 0 n Pl .l/3 Gj .1/4 

(l + 0. 494 Pr 2 / 3 ) 1 ^ 4 

for a value of Pr between 0.74 and 0. 90. Equation (B7) is similar in form to the empir- 


(B7) 
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ical relation for free convection flow along vertical planes and cylinders given in refer- 
ence 12 as 


Nu 


- C,ljL - n(GrPr) m 
k 


(B8) 


Equation (B8) is used herein with exponent m = 1/3 and a multiplier n = 0. 14 because it 
is simpler and fits the data somewhat better than equation (B7). 

At this point, some choice of Tg that is consistent with the definition of h c ^ and 
fits the data for q^g (i*e. > ^ (U^)) must be made. In reference 13, liquid hydrogen 
was pressurized with a low mixing diffuser and no liquid outflow. The adiabatic com- 
pression temperature given by T a( j = T q (P/P o )^”^^ y was used as the choice for Tg. 
This relation gave good agreement between analytical and experimental mass flux results 
for hydrogen pressurant. For the conditions described in reference 13, appreciable con- 
densation occurred. Reference 13 indicates that as Tg increased, there was a tendency 
toward evaporation, away from the condensation results that occurred when the adiabatic 
temperature was used. With greater ullage gas mixing (due to diffuser characteristics as 
well as the liquid outflow process), Tg would be expected to be greater than T af j. For 
the work described in this report, Tg was obtained from the following relation evaluated 
with experimental data: 

^U-S = h c, L (T 5 ‘ T sat ) = “ U L 


One experimental run, where the expulsion time was 271 seconds (run 95), was used 
for the determination of Tg where h £ ^ is given by (B8). For this condition Tg was 
determined to be 1. 3 times the adiabatic compression temperature or 46.3 K (83. 3° R) 
for a tank pressure of 34. 47x10^ newtons per square meter (50 psia). For convenience 
this value of Tg was used for all comparisons since the effect of the interfacial terms in 
the energy equation was small in the experimental situation. 

Using equation (B8), as well as the values for T g and T gat discussed previously, 
gave the final form of the equation used to evaluate the heat transferred from the ullage < 

gas to the liquid interface as follows: 

Op— S = - (°- 14)(GrPr) 1/3 (46. 3 K - T gat ) (B9) 

L 


In order to incorporate liquid heating to the analysis the term 
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A L _ k(0 - 14)(GrPr) 1/3 (46. 3 K - T gat ) • A L 

X^p LX^Vp 

must be added to the right side of equation (Al). 


(BIO) 
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APPENDIX C 
RAMP ANALYSIS 

The amount of pressurant gas needed to pressurize a propellant tank initially may be 
important for multiburn missions. When the coast period between firings is long, a col- 
lapse in ullage pressure may develop. Under these circumstances a significant amount of 
pressurant may be necessary to repressurize the propellant tank for the next firing. 

An application of the work reported in reference 1 is the prediction of the pressurant 
requirements for the initial pressurization (ramp) and hold periods. A separate com- 
puter program, which determines the pressurant as well as energy requirements during 
the ramp and hold periods, is described herein. The same equations which describe the 
expulsion period, as outlined in appendix A, are also applicable for the ramp and hold 
periods. However, the ramp period is more difficult to model analytically than the ex- 
pulsion period, particularly for small ullage volumes where the liquid surface is near the 
injector outlet. Even though the experimental results indicated relatively large amounts 
of mass transfer during this period, the incorporation of mass transfer into the ramp 
analysis was not attempted because of the added complexity and incomplete knowledge of 
the mass transfer phenomenon. The heat transfer from the ullage to the liquid surface is 
also neglected in the analytical model for the ramp period. 

The analysis used for the ramp period computes the gas temperatures in the ullage at 
any time during the pressure rise from the gas energy equation. The corresponding gas 
velocities are computed from the equation of continuity. The iterative method to be de- 
scribed shows how convergence is achieved in the solution of the gas energy and continu- 
ity equations. The predicted mass of pressurant is based on an integration of the volume 
elements in the ullage at the end of the ramp and hold periods assuming a one component 
ullage (100 percent helium). Quantitatively, the entire mass of pressurant requirements 
for the ramp period was less than the expulsion period by a factor of 30 to 60 when the 
initial ullage was 4 percent of the tank volume. 


INPUT DATA REQUIREMENTS 

For the solution to proceed, a set of boundary and initial conditions are required. 
These conditions, which are the same for the expulsion as well as the pressurization, are 
as follows: 

(1) At time t = 0, the values of gas temperature T and wall temperature T w as 

functions of x, the position within the ullage 

(2) On the boundary x = 0, the value of inlet gas temperature T as a function of time 
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(3) At the liquid surface, the value of gas temperature T, wall temperature T , and 
velocity V as functions of time (Although movement of the interface has been 
noted during the ramp pressurization period, no significant effect on the pro- 
gramed output was noted with the value of V = 0 at the interface. ) 

(4) Tank pressure P, outside heating rate q w , inside hardware heating rate q^ as 
a function of time (Like the other initial conditions, the pressure P as a function 
of time or ramp pressure curve is defined by a discrete set of points which ap- 
proximate a smooth curve. In regions of pronounced curvature, more points are 
needed for accurate definition than for linear portions. ) 

(5) Constant value of heat-transfer coefficient h c , or a correlating equation from 
which h c may be evaluated at each net point from values of T, T w , and P 

(6) Tank radius as function of axial distance down from the top of the tank 

(7) Tank wall material properties: density p,„ and specific heat C,„(T ) 

(8) Tank wall thickness (average membrane plus weld area thickness) as function of 
axial distance down from the top of the tank 

(9) Pressurizing gas properties: Molecular weight M, specific heat C (T), and 

ir 

compressibility factor Z(P, T) 

(10) Initial ullage height, total time of run, the number of net points in the initial ul- 
lage space 

(11) The initial time step At used in following the pressure rise as well as establish- 
ing the points of computation 

(12) If the hold period is to be included in the analysis, then the time for the end of the 
ramp must be specified 


APPLICATION OF BASIC EQUATIONS 

Reference 1 makes the substitution of T w ^ from the finite difference form of equa- 
tion (A6) into the finite difference form of the energy equation. Rearranging gives a 
quadratic in the gas temperature T) where the prime refers to a step forward in time 
and the quantities without the prime are evaluated at the previous time step: 



The quantity marked with the asterisk may be evaluated either at the beginning or the end 
of the time interval. 

A difficulty can arise when evaluating the gas energy equation expressed by the pre- 
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vious quadratic. This occurs when the heat transfer takes place from the wall into the 
ullage gas. For this situation, the solution of the continuity equation provided negative 
gas velocities which made it impossible for equation (Cl) to converge on the real roots. 

At the start of the ramp (immediately after filling the tank), the initial wall tempera- 
ture distribution in the ullage is higher than the gas temperature distribution. This is 
brought about since the wall surface above the liquid is exposed to the ambient tempera- 
ture. But the ullage gas temperature near the liquid interface is close to the saturation 
temperature at 1 atmosphere. 

The technique used when T . > T^ involved a direct substitution. The finite dif- 
ference form of equation (A4) is 



Combining the equations involving ullage gas temperature T? and velocity V?, equa- 
tions (Cl) and (C2), gives the following cubic results: 


Vi + 


{tI + < t {t\ % T * ■ < t T ‘-i ^ + < T * t b *_ 
-<£(Ti% T rt-'-< T $\ T ™- e 


T? 

1 


(C3) 


This arrangement of terms eliminates V| and the cubic is solved for the gas tem- 
perature T|. 


ANALYTICAL PROCEDURE 

The analytical procedure uses a variable time increment At in following the pres- 
sure rise. With this technique, the iteration was stable over a range of inlet conditions 
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and the results were consistent with the recorded data. 

First, the velocity distribution is determined in the initial ullage at time t = 0. The 
substituted form of equation (A4) (in which the energy equation (Al) is substituted for 
3T/3t) is put into finite difference form and solved for the velocity distribution for the 
first time. 

For most of the ramp rims encountered in this investigation, an initial time incre- 
ment of 1 second proved satisfactory. 


TEMPERATURE CALCULATIONS FROM TOP TO INTERFACE 


Since values of V have been obtained at each net point at time t = t^ = 0, attention 
is turned to equation (C3) which is cubic in T?. During the iteration, the cubic equa- 
tion (C3) is solved for the gas temperature T? starting at the point N 2 in figure 21. 


Inlet gas M, T(inlet) at Nj 



Figure 21. - Analytical model. (Coordinate system is positive in downward direction from 
x=0 at Nj to x=n at interface N z .) 
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When this equation is first solved for the ullage temperature distribution, a value for 
T? + j is not available- A substituted value of T i + 1 proved to be satisfactory as an ini- 
tial guess to get convergence. All other quantities in equation (C3) are available from 
the initial conditions. The value for T?_j, the temperature at Np is known as a bound- 
ary condition. 

The solution for T? at Ng follows, and this procedure continues to calculate gas 
temperatures until the boundary at the interface is reached. The values for the corre- 
sponding wall temperatures are calculated using the finite difference form of equa- 
tion (A6). 


VELOCITY CALCULATIONS FROM INTERFACE TO TOP 

Although the ullage temperatures are computed starting at the top (fig. 21), the ve- 
locity equation (C2) is used to calculate the ullage gas velocity starting with the point N y 
near the interface. The velocity at the interface N , the boundary value, is zero with 
no expulsion. 

The ullage gas velocity is calculated from point to point until the top of the tank is 
reached. The new velocities are used in equation (C3) along with previous values of 
t! + i and the temperature distribution is redetermined. This process is continued until 
convergence is achieved over the entire ullage. The time is advanced to tg and a new 
set of velocities is determined. 

COMPLETING THE SOLUTION 

With the new velocities at time tg, we evaluate equation (C3) again starting at point 
Ng and terminating at the interface. A value for Tj + j is always available from the 
previous iteration, although a substituted value of T. + j is used as the first value. 

The new values for T! at all the points for time tg are used to recompute the ve- 
locity distribution. This new set of velocities is then compared with the previous set and 
convergence is assumed if the deviation is less than half of 1 percent for every velocity 
in the time set. A time step is then taken to tg. 

If convergence is not achieved after 40 iterations, the time step is reduced and the 
iteration process is reinitiated. Generally the reduction in time step becomes necessary 
only when there is a severe change in the slope of the ramp curve particularly in the 
early stages of the pressure rise. 

For the new time tg, the temperature T! in equation (C3) is determined from the 
converged value using the iterating method. This procedure continues to evaluate the 
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gas temperatures and velocity distribution in the ullage for each time step taken in fol- 
lowing the rate of pressure in the tank. 

The initial gas velocity distribution used in solving equation (C3) for each new time 
t is obtained from the previous time as follows: 



This iterative procedure can be used for a constant pressure representing the hold 
period. However, for initiating the ramp, an actual pressure rise must be used. A typ- 
ical example is shown in figure 22. 



Time, sec 

Figure 22. - Tank pressure as function of time during 
initial pressurization period for run 96. 
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APPENDIX D 


DATA REDUCTION 


PHYSICAL DESCRIPTION OF PROBLEM 


An initially vented tank containing two-phase hydrogen was pressurized from 1 at- 
mosphere to a new pressure by adding helium gas. The system was then allowed to sta- 
bilize at which time liquid outflow was started. During this expulsion period, pressurant 
gas (helium at constant temperature) was added to the tank at a rate that maintained a 
constant tank pressure while expelling the liquid at a desired rate. The amount of pres- 
surant gas used during the expulsion phase is dependent on (1) the volume of liquid dis- 
placed with no heat or mass transfer, (2) the heat transfer to the tank wall and liquid, 

(3) the amount of mass condensed or evaporated, and (4) the amount of helium gas ab- 
sorbed by the liquid hydrogen. 

The main parameter used in the comparisons was the ratio of the ideal pressurant 
requirement to the actual pressurant requirement. The ideal pressurant was determined 
under the assumption that the incoming pressurant gas did not exchange energy or mass 
with the surroundings. Under this assumption, the ideal pressurant required for the ini- 
tial pressurization of the tank was determined by the adiabatic relation 


M I = 


MP V 
o o 

ZRT~ 



(Dl) 


The ideal pressurant required for the expulsion period was determined by the relation 

MP AV tt 

M t = y. (D2) 

1 ZRT g 


MASS BALANCE 


A mass balance was performed on the ullage volume from an initial time t^ to a 
final time tj as follows: 


M U,f " M U,i + 


M G, i— f ± M t, i— f 


(D3) 


A discussion of how the terms of equation (D3) were determined follows. 
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Pressurant Gas Added (Mq j_^) 


The weight of the actual pressurant gas added from any initial time t to any final 
time tj was determined by numerical integration of the gas orifice equation 

m G = f tf YD 2 C Vp AP* dt (D4) 

A 

The time necessary to expel AVy of liquid is (tj - tj). 


Ullage Mass 

The initial ullage mass My . and final ullage mass My ^ were obtained by numer- 
ical integration of the particular density profiles: 


M U, i ■= A. . P dv “ E Vn where P - f ( T ’ X H,’ X He) <° 5 > 

U>1 n=l \ i / 

Nf . » 

M U, f = fv v f P dV “ E P n V n Where P = f ( T ’ X H 2 ’ X He) (D6 > 


The internal tank volume was considered as 57 (corresponding to thermopile loca- 
tion) horizontal disk segments. Each of these segments was in turn divided radially into 
a series of concentric rings, the number of which depended on location of radial temper- 
ature sensors and the vertical position of the disk segment being considered. These 
rings (339 in all) comprised the V n 's in the previous calculations. In this manner, 
vertical temperatures as well as radial temperature gradients could be incorporated into 
the mass calculations. The position of the liquid level prior to and after expulsion de- 
termined the number of gas volume rings (N^ and N^) used in the ullage mass calcula- 
tions. The density of the two component ullage was calculated from the Beattie- 
Bridgeman equation of state for a mixture given in reference 14 as 

p = RT(1 -e _ ) (y + B ) _ (D7) 

ITT* 2 1 2 


where 
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l T " ( X H 2 V A H 2 + X He V A He) 2 


X H 2 a H 2 + X He a He 


) 


4 


B T - \ X H 2 B H 2 + X He B He / 11 


X H 2 b H 2 + X He b He ' 


and 


e = • 


X H 2 C H 2 + X He C He 
vT 3 


Equation (D7) can then be solved for ^ (where 4- = p n ) by knowing the pressure, tem 


n n 


perature, and ullage gas concentration for each AV n segment. 


Mass Transfer 


The mass transfer was calculated from equation (D3) as a result of knowing £, 
My j and M G that is, 


M t, i—f _ M U, i + M G, i— f " M U, f 


(D8) 


If M t . ^ was a positive quantity, mass was considered leaving the ullage volume (e.g. , 
condensation and/or absorption). 


ENERGY BALANCE 

For the thermodynamic system consisting of the entire tank and its contents (tank 
+ ullage gas + liquid), the first law of thermodynamics for an increment of time dt may 
be written as 


dU T - (6 M g )(u g 


P G V G 


H 

2g 


+ Z 


Cl ~ ( 5M L ) (u L + P L V L + + z L j + 5Q - 5W 


(D9) 


The kinetic and potential energy terms are small in comparison with the other en- 
ergy terms and are neglected in this development. If h = u + Pv is substituted, equa- 
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tion (D9) becomes 


dU T = (6M G )h G - (6M L )h L + 6Q - SW (DIO) 

For this system, there is no external work done so 6W = 0 and the final form of equa- 
tion (D9) therefore becomes 


dU T = (SM G )h Q - ($M L )h L + 6Q 


(Dll) 


Equation (Dll) can be integrated over any time period. The physical interpretation 
of the quantities in equation (Dll) is as follows: 



Change in 
system 
energy 
(Tank + Gas 
+ Liquid) 


£ 


f M G h G dt 


Energy input 
by pressurant 
gas inflow 


£ S M L h L dt 

^ — V ' 

Energy leaving 
through 
liquid outflow 



Energy from 
environment 


(heat leak from 
conduction, 
convection and 
radiation) 


(D12) 


A discussion of how the terms of equation (D12) were evaluated follows. 


Energy Input by Pressurant Gas in Flow 


The first term in equation (D12) may be evaluated as follows: 


n=- 


tf-ti 

At 


J f M G h G dt ~ ^2 ^G, n h G, n At 

l n-0 


(D13) 


The pressurant flow rate M G was determined from equation (D4). The specific en- 
thalpy of the inlet gas was evaluated at the inlet temperature and pressure at each time 
increment At. 
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Energy Leaving by Liquid Outflow 


The energy of the liquid that leaves the system can be evaluated as follows: 



M L h L 


dt — 


n- 


V^i 

At 


E M t h T „ At 
L, n Li, n 

n=0 


(D14) 


The liquid flow rate was determined from the turbine flowmeter. The specific en- 

thalpy of the liquid was evaluated at the outlet temperature. 


Energy Input from Environment 

The rate of energy input into the tank from the environment was assumed to be the 
same for all cases and was determined from a boiloff test. This test indicated a nominal 
value of 0. 732xl0 3 joules per second (0. 694 Btu/sec) should be used. This value, which 
includes heat input by radiation, convection, and conduction through pipes and supports, 
was in all test cases less than 7. 0 percent of the energy added to the tank by the pres- 
surant gas: 



Q dt « 0. 732X10 3 


(t f 


l i) 


(D15) 


Change in System Energy 

The change in system energy can be separated into three categories: (1) change in 
ullage energy, (2) change in liquid energy, and (3) change in the wall energy. Stated 
mathematically, 

dU T = dUy + dU L + dU w (D16) 

Change in Ullage Energy 

The change in the ullage energy over any given time interval (T — tj) is obtained by 
subtracting the internal energy of the ullage at time t from the internal energy at time 
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Making use of the relation U = H - PV gives 




(D18) 


The ullage gas density was determined using equation (D7). The ullage gas enthalpy 
was calculated using the relation 


h Mh 2 k M He h 

h — — — IItt h h-rj _ 

M t n 2 M t e 


(D19) 


Therefore, since the pressure and temperature and ullage gas concentration profiles at 
times tj and t^ were known, the change in ullage energy was evaluated. 


Change in Liquid Energy 


The change in energy of the liquid in the tank can be determined in a manner similar 
to the change in ullage energy: 



where the liquid density and enthalpy are functions of pressure and temperature. 



Change in Wall Energy 


The change in wall energy was determined by applying the first law of thermody- 
namics to an element of the wall: 


r\ 

\ 


f 


dU w =AU w = AM w # C y dT where C y = C y (T) (D22) 


The total change of the wall is then 


iU w,T S E 

M w 


AU = 
w 


T 

' E iM w L 2 c v< T > dT 

M w 1 


(D23) 


Total Energy Change of System 


For convenience equation (D16) is substituted into (D12): 

/ tf i [UlJ + “w + UL]d ‘ ' j( t£ “ GhG * " _[ i{ “L h L * + ^ « 


dt (D24) 


Rearranging terms gives 


£* [M G h G * Ql® = £* I“L h L dt + dU lJ + C idV U + jC* 

> -i .. > 


dU w (D25) 


Total energy Total change in liquid Total change Total change 
added in tank plus liquid in ullage in wall energy 

(AUrp) expelled energy energy 

(AU l ) (AU n ) 


< AU w> 


Dividing through by AU,p gives 


± _ AU L , AU U , AU w 
AU t AU t AU t 


(D26) 
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The data presented herein are in the form of these ratios which show the relative distri- 
bution of the total energy input. 
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